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0.  Background 


This  report  is  one  of  a  series  of  studies  stemming  from 
current  research  on  electrohydrodynamic  (EHD)  power  genera¬ 
tion,  sponsored  by  the  Department  of  Energy.  A  certain  EHD 
scheme  proposed  by  the  Marks  Polarized  Corporation  is  presently 
under  development.  The  authors  of  this  report  have  partici¬ 
pated  as  consultants  on  various  theoretical  and  analytical 
aspects  of  this  effort.  One  of  our  main  objectives  has  been 
to  make  an  independent  appraisal  of  the  performance  possi¬ 
bilities  and  limitations  inherent  in  the  general  type  of  EHD 
system  involved  in  this  development.  Another  objective  has 
been  to  create  a  sound  analytical  basis  for  optimizing  the  key 
design  parameters  of  such  a  device,  and  to  suggest  design 
improvements  where  appropriate. 

The  results  of  our  first  study  are  reported  in  detail  in 
Ref.  (1)  and  are  summarized  more  concisely  in  Ref.  (2) .  This 
work  indicated  that  the  design  as  initially  proposed  would 
yield  unacceptably  poor  performance  and  suggested  a  certain 
design  revision.  Further  work  was  done  to  analyze  such  a 
revised  design  and  is  reported  in  the  appendix  of  Ref.  (1) . 

An  improved,  simplified  and  shortened  version  of  this  work  is 
reported  in  Ref.  (3).  Unfortunately,  our  results  indicated 
that  the  revised  design,  while  showing  somewhat  improved  per¬ 
formance,  still  would  not  attain  competitive  overall  efficiency. 

The  Marks  Polarized  Corporation  has  challenged  these 
pessimistic  conclusions.  They  argue  that  in  our  analysis 
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the  compressibility  of  the  medium  was  not  taken  into  account 
and  that  if  it  had  been,  the  predicted  performance  would  have 
been  drastically  improved.  We  have  good  reasons  to  conclude 
that  Marks'  argument  about  compressibility  is  fallacious.  In 
the  first  place  our  analysis  did  consider  all  of  the  major 
effects  of  compressibility;  only  certain  minor  effects  in  one 
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part  of  the  system  were  neglected  in  order  to  simplify  the 
calculations.  While  a  more  elaborate  and  accurate  treatment 
of  these  neglected  effects  is  certainly  possible,  we  maintain 
that  this  refinement  will  not  suffice  to  change  the  generally 
pessimistic  conclusions  established  in  our  original  analysis. 

Fortunately,  this  difference  of  opinion  between  ourselves 
and  the  Marks  Polarized  Corporation  pertains  to  a  question  of 
'fact  anout  which  it  is  entirely  possible  to  obtain  objective 
evidence.  To  this  end  we  present  in  this  report  a  revised 
and  improved  analysis  which  takes  all  compressibility  effects 
exhaustively  and  minutely  into  account.  Calculations  based  on 
this  revised  analysis  should  eventually  settle  the  above 
question  one  way  or  the  other.  A  typical  trial  calculation  is 
presented  later  in  this  report.  This  initial  calculation  has 
been  made  manually,  but  systematic  further  work  will  be  done  by 
computer.  The  full  results  and  conclusions  derived  from  these 
calculations  will  be  reported  later. 

It  has  also  been  conjectured  by  Marks  and  others  that  it 
might  be  possible  to  improve  greatly  the  performance  of  the 
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EHD  generator  by  employing  a  two-fluid  cycle.  We  therefore 
present  in  this  report  a  detailed  analysis  and  procedure  by 
means  of  which  the  performance  of  any  specified  two-fluid  EHD 
system  can  be  calculated.  Our  analysis  takes  into  account 
fully  the  compressibility  of  both  fluids  at  all  locations  in 
the  system. 

We  propose  to  undertake  systemic  computer  calculations 
of  this  kind  which  should  enable  sin  objective  assessment  to 
be  made  of  the  performance  potentialities  and  limitations  of 
several  proposed  two-fluid  EHD  generators.  Such  results  would 
be  of  great  value  in  helping  decide  whether  the  two-fluid  EHD 
generator  shows  sufficient  promise  to  warrant  further  develop¬ 
ment  effort  and,  if  so,  what  the  optimum  design  point  values 
of  the  key  parameters  should  be. 

As  far  as  we  know,  the  present  report  is  the  only  pub¬ 
lished  analysis  of  the  two-fluid  EHD  generator  which  includes 
the  compressibility  of  both  fluids  and  which  considers  the 
complete  thermodynamic  cycle.  The  principal  earlier  effort 
along  these  lines  is  given  in  Ref.  (4) ,  but  the  present  report 
goes  well  beyond  anything  undertaken  in  Ref.  (4) .  In  particu¬ 
lar,  Ref.  (4)  makes  no  attempt  to  analyze  the  complete  thermo¬ 
dynamic  cycle.  Moreover,  although  it  contains  certain  useful 
data  and  is  correct  in  some  respects.  Ref.  (4)  is  flawed  by  a 
number  of  serious  errors. 
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While  the  present  report  is  otherwise  self-contained,  it 
does  assume  that  the  reader  has  a  general  familiarity  with  the 
type  of  EHD  power  generator  described  in  Refs.  (1),  (2)  and 
(3) .  It  is  recommended  that  the  reader  unfamiliar  with  this 
general  scheme  first  review  briefly  the  introductory  material 
in  the  first  section  of  any  one  of  these  references  before 
proceeding  with  the  more  detailed  treatment  in  this  report  as 
this  earlier  background  information  is  not  covered  or  repeated 
in  the  present  report.  We  do,  however,  include  here  a  schema¬ 
tic  diagram.  Fig.  0.1,  which  is  taken  from  Ref.  (3).  It 
illustrates  one  particular  version  of  an  EHD  power  generator. 


PUMP 


1. 


Introduction 


For  purposes  of  analysis,  the  overall  thermodynamic  pro¬ 
cess  in  the  two-fluid  EHD  generator  may  be  represented  by  the 
schematic  flow  diagram  shown  in  Fig.  1.1. 

Primary  fluid  enters  the  ejector  at  station  1,  secondary 
fluid  enters  at  station  2  and  the  mixture  leaves  at  station  3. 

The  mixture  then  passes  through  a  diffuser,  process  3-*-4, 
the  basic  purpose  of  which  is  to  reduce  the  kinetic  energy  of 
the  flow  and  thereby  decrease  the  associated  friction  losses 
downstream  of  the  diffuser.  On  the  other  hand  the  diffuser 
introduces  certain  losses  of  its  own  so  that  some  care  is 
necessary  to  optimize  the  design  in  this  respect. 

The  gas  mixture  then  passes  through  the  electrical  power 
section,  process  4 -*-5 ,  from  which  the  gross  electrical  power  is 
extracted. 

Finally,  the  gas  mixture  enters  the  condenser/separator 
at  station  5.  The  primary  fluid  is  separated  from  the  gas 
mixture  by  condensation  and  leaves  at  station  6.  The  fluid  at 
this  point  is  assumed  to  be  compressed  liquid  at  a  known  static 
pressure  Pg  and  a  known  static  temperature  Tg  .  Kinetic 
energy  at  station  6  is  regarded  as  negligible.  The  secondary 
fluid  leaves  the  condenser/separator  and  enters  the  ejector  at 
station  2.  We  assume,  moreover,  that 

P6  =  P2  =  Px  (1.1) 

and 

T6  =  T2  (1.2) 
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Fig.  1.1  SCHEMATIC  FLOW  DIAGRAM  OF  TWO -FLUID  EHD  CYCLE. 


The  primary  fluid  which  leaves  the  condenser/separator 
at  station  6  is  circulated  back  through  the  pump  and  through 
the  boiler/superheater,  process  7+1. 

It  should  also  be  explained  that  the  primary  fluid  at 
station  1  is  a  condensible  vapor  of  high  molecular  weight 
which  contains  fine  liquid  droplets.  Unless  specifically 
stated  otherwise  the  quality  z  at  this  location  is  taken 
as  0.95.  The  important  thermodynamic  static  properties  at 
station  1  can  be  expressed  in  the  following  form. 

h^  =  h^  +  z  h^g  Enthalpy  (1.3) 

s1  =  sf  +■  z  Sf  Entropy  (1.4) 

Here  subscript  f  denotes  the  saturated  liquid  and  subscript 
fg  denotes  the  change  between  saturated  liquid  and  saturated 
vapor.  All  properties  are  evaluated  at  a  specified  pressure 
;  the  corresponding  saturation  temperature  T^  is  then 
also  known. 


The  corresponding  stagnation  properties  Tg^  and  hg^ 
are  fixed  when  the  stagnation  pressure  Ps^  is  specified. 

Of  course  the  entropy  has  the  same  value  s^  at  the  stagna¬ 
tion  condition  as  at  the  corresponding  static  state.  Thus 
specification  of  Pg^  and  s-^  suffices  to  fix  all  other 
stagnation  properties. 

While  it  is  convenient  to  start  the  analyses  by  stipula¬ 
ting  the  quality  at  station  1,  it  is  also  advisable  later  in 
the  calculation  to  determine  the  corresponding  quality  at 


station  4,  the  entrance  to  the  electrical  working  section. 

It  is  essential  that  the  fluid  entering  the  working  section 
contain  finely  dispersed  liquid  droplets  of  the  proper  size. 
These  droplets  carry  the  electrical  charges  that  are  the 
central  feature  of  EHD  power  generation.  We  assume  tentatively 
that  acceptable  droplet  characteristics  are  obtained  by 
maintaining  the  quality  at  the  inlet  to  the  working  section 
in  the  range  0.92  to  0.98.  These  limits  are  only  estimates. 

The  thermodynamic  analysis  would  be  greatly  simplified 
if  both  fluids  were  perfect  gases  with  constant  specific 
heats.  In  fact  only  the  secondary  fluid  satisfies  this  re¬ 
quirement.  Nevertheless,  we  can  define  an  "equivalent  perfect 
gas"  which  adequately  approximates  the  essential  thermodynamic 
properties  of  the  primary  vapor/liquid  mixture.  When  carried 
out  judiciously,  this  procedure  simplifies  the  analysis  at  the 
cost  of  only  a  small  loss  in  numerical  accuracy^. 

The  ratio  of  specific  heats  y  of  the  "equivalent  perfect 
gas"  may  be  deduced  from  the  important  relation 

Yl/(Yl  -  1) 

<p3l/pi>  *  <Tsl/Tl>  d.5) 


Solving  this  for  y^  gives 


Y1  = 


1  - 


in(Tgl/T1) 


,-l 


(1.6) 


Another  fundamental  property  of  the  "equivalent  perfect 
gas"  is  its  specific  heat  Cpl  which  may  be  evaluated  from 
the  expression 


(1.7) 


'Pi 


<hsl  -  hl> 
(Tsl  -  Tl> 


where  all  quantities  on  the  right  are  now  known. 

It  then  follows  from  standard  perfect  gas  relations  that 
the  gas  constant  and  the  molecular  weight  of  the 

"equivalent  perfect  gas"  are,  respectively, 


R1  " 


'*1  "  ^ 


Si 


(1.8) 


and 


where  R  =  universal  gas  constant 
=  8315  joules/kg  °K 

Eqs .  (1.6)  through  (1.9)  complete  the  definition  of  the 
"equivalent  perfect  gas"  which  adequately  simulates  the  essen¬ 
tial  thermodynamic  properties  of  the  primary  vapor/liquid 
mixture. 

It  should  also  be  explained  that  the  secondary  fluid  at 
station  2  is  predominantly  a  noncondensing  gas  of  low  molecu¬ 
lar  weight  primary  vapor. 

It  can  be  shown  that  the  mass  ratio  oi  of  condensible 
primary  vapor  to  dry  secondary  gas  at  station  2  is  given  by 
the  expression 


01 


W 


WG  ^ 


w 
-  vv J 


(1.10) 
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where 


Wv  =  true  molecular  weight  of  condensible  vapor  (not  to 
be  confused  with  the  "equivalent  molecular  weight 
W^"  considered  earlier) 

WG  =  molecular  weight  of  noncondensing  gas 
Pv(Tg)  =  vapor  pressure  of  condensible  primary  fluid  at  known 
temperature  Tg 

0 

P2  =  known  static  pressure  at  station  2 

The  present  analysis  is  restricted  to  conditions  under 
which  a)  is  very  small  compared  with  unity  and  may  be  neglec¬ 
ted.  Thus  the  fluid  at  station  2  may  be  treated  as  dry  gas. 
This  greatly  simplifies  the  analysis.  Eq.  (1.10)  is  useful  for 
verifying  that  this  assumption  is  indeed  satisfied  in  any  par¬ 
ticular  instance. 

It  should  be  added  that  if  it  later  becomes  of  interest 
to  make  calculations  for  circumstances  under  which  parameter 
to  is  not  negligible,  the  present  analysis  can  in  fact  be 
generalized  to  include  this  effect.  For  the  present,  however, 
we  prefer  to  deal  only  with  the  simpler  situation  in  which  o> 
may  be  neglected. 


2.  Ejector 

The  ejector  receives  a  primary  stream  of  high  molecular 
weight  gas  at  station  1,  a  secondary  stream  of  low  molecular 
weight  gas  at  station  2  and  discharges  the  resulting  mixture 
at  station  3.  Static  conditions  at  these  three  stations  are 
designated  by  subscripts  1  ,  2  ,  3  .  Stagnation  conditions 

at  the  corresponding  stations  are  designated  by  subscripts 
si  /  s2  and  s3  . 

The  following  quantities  are  arbitrarily  specified  or 
known  at  stations  1  and/or  2,  namely, 

Molecular  weights:  and  W2 

Ratios  of  specific  heats:  and  Y2 

Stagnation  pressure  of  primary  jet:  P  , 

S  X 

Stagnation  temperature  of  primary  jet:  T  , 

Static  pressures:  Pjl  *  P2  (=  Pg) 

Static  temperature  of  primary  jet:  T^  (=  saturation  tempera¬ 
ture  of  primary  fluid  at  pressure  P^) 

Static  temperature  of  secondary  jet:  T2  (=  Tg) 

Mass  flow  ratio:  =  m^/Cn^  +  ih2)  =  x 

Velocity  ratio:  V]/V2  =  ^ 

We  treat  the  fluids  at  stations  1,  2,  3,  4  and  5  as  per¬ 
fect  gases  with  constant  specific  heats. 

The  velocity  ratio  of  the  ejector  may  be  developed  as 
follows: 


i  Mi  Ri  Tsi  ,  y  .  _i 

a2  M2  y  Y2  R2  T2  f77h  M2 


(2.1) 
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We  can  now  use  this  result  and  other  standard  perfect 
gas  relations  to  complete  the  calculation  of  key  thermody¬ 
namic  properties  at  stations  1  and  2.  It  is  convenient  to 
arrange  these  calculations  sequentially  as  follows : 


M 


1 


M 


2 


Ts2  *  T2  (Ts2/T2) 


(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 

(2.7) 


Before  we  can  determine  the  actual  conditions  at  station 
3,  the  outlet  of  the  ejector,  it  is  first  necessary  to  analyze 
two  other  hypothetical  cases  as  indicated  schematically  in 
Fig.  2.1.  In  all  three  of  these  cases  the  flow  is  taken  as 
steady  and  adiabatic. 

The  first  case,  that  shown  in  Fig.  2.1(a) ,  involves  a 
hypothetical  device  which  receives  two  gas  streams  at  stations 
1  and  2  and  discharges  the  resulting  mixture  at  station  x  .  In^ 
let  conditions  at  1  and  2  are  identical  to  those  of  the  actual 
ejector. 
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Recall  that  the  adiabatic  mixing  of  two  different  gases  is 
an  inherently  irreversible  process  that  always  involves  a  cor¬ 
responding  entropy  increase.  We  stipulate  that  the  first  case 
is  such  that  the  only  irreversibility  which  occurs  is  that  asso¬ 
ciated  with  this  mixing.  We  term  this  mass  mixing  to  distin¬ 
guish  it  from  another  type  of  mixing  considered  below  which  we 
term  momentum  mixing.  * 

The  second  case,  that  shown  in  Pig.  2.1(b),  involves  an 
ideal  ejector  which  receives  two  gas  streams  at  stations  1  and 
2  and  discharges  the  resulting  mixture  at  station  y.  Again, 
the  inlet  conditions  at  1  and  2  are  identical  to  those  of  the 
actual  ejector.  This  case  satisfies,  among  other  relations,  the 
idealized  one-dimensional  momentum  equation  for  frictionless, 
constant  area  flow.  It  can  be  shown  that  in  this  case  there  is 
an  overall  entropy  increase  which  includes  not  only  the  previ¬ 
ously  mentioned  effect  of  mass  mixing  but  also  a  further  increase 
associated  with  momentum  mixing. 

The  third  case,  that  shown  in  Pig.  2.1(c) ,  represents  the 
actual  ejector  itself  which  receives  input  streams  at  stations 
1  arid  2  which  discharges  the  resulting  mixture  at  station  3. 

All  three  of  the  above  cases  refer  to  constant  area  pro¬ 
cesses  in  the  sense  that 

Ax  =  Ay  =  A3  =  (Ax  +  A2)  (2.8) 

Once  the  hypothetical  operating  conditions  at  stations  x 
and  y  have  been  found,  it  then  becomes  possible  to  express  the 
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corresponding  actual  conditions  at  station  3.  This  is  accom¬ 


plished  through  the  use  of  an  ejector  effectiveness  nE  . 

This  parameter  is  defined  more  explicitly  in  the  later  analysis. 
Values  of  nE  must  be  estimated  by  reference  to  test  data  on 
ejectors  (i.e.,  Ref.  6,  7,  8). 

9 
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3.  Adiabatic  Mixing  of  Two  Streams 


The  process  in  the  ejector  is  treated  as  the  adiabatic 
mixing  of  two  streams.  The  following  relations  can  be  shown 
to  apply. 


R  -  8315  Joule/kg  K 


R_  = 

"l  W1 


R  -  R 
2  ~ 


R„  =  R  =  R,  =  X  R.  +  (1  -  x)  R„ 
X  y  3  1  ^ 


Yi 


'pi  (y1  -  1) 


Y2  R2 

'P2  =  (Y,  -  1) 


:  =  c  =  c  ,  =  x  c  .  +  (l-x)  c  _ 
px  py  p3  pi  p2 


IeI 


Yx  =  Yy  =  Y3  =  icp3  -  R3> 


rsx=Tsy  =  Ts3=  Ix  Cpl  Tsl+  (1 


Ax  =  Ay  =  A3  =  (A1  +  V 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 


(3.6) 

(3.7) 


x)  Cp2  Ts2] 


(3.8) 

(3.9) 

(3.10) 


Notice  that  the  quantities  R  ,  C  ,  y  r  T  and  A  all 

p  s 

have  the  same  values  at  stations  x,  y  and  3. 

This  fact  may  be  used  to  simplify  many  of  the  subsequent 
expressions.  It  is  always  permissible  to  substitute  subscript 
3  for  subscripts  x  or  y  on  any  of  these  quantities. 
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4.  Mass  Flow 
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also 


mx  = 


Psx  A3 
VR 3  Ts3^3 


f(Mx) 


K  =  — sy . 3  f(M 


"VR3  Ts3/y3 


*3  - 


Ps  3  A3 
VR3  Ts3/y3 


f  (m3) 


'*l\ 

ST )  =  x 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


*x  =  *y  =  *3  =  *]_  +  *2  (4.10) 

From  Eqs.  (4.4),  (4.5),  (4.9)  and  (4.10)  we  readily  find 

that 


/Ps2\  /R1  Tsl  Y2  £(M2} 
\Psl/v  R2  Ts2  y1  f  (M1) 


(4.11) 


Then  from  Eq.  (2.10)  we  infer  also  that 


Ax\  (Aj/Aj) 

A 3/  =  (1  +  A1/A2) 


(4.12) 


Since  all  quantities  on  the  right  side  of  Eq.  (4.11)  are 
now  known,  Eqs.  (4.11)  and  (4.12)  fix  the  area  ratios  Aj/A2 
and  aj/a3  • 

From  Eqs.  (4.5),  (4.6),  (4.9),  (4.10)  and  (2.8)  we  deduce 
further  that 
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....  ,  _  (1  -  W  /PS2 

£(Mx»  -  n-^r-  fc. 


/R3  Ts3  Y2 
3  S3  f  (M  ) 

2  s2  y3  Z 


(4.13) 


All  quantities  on  the  right  side  of  Eq.  (4.13)  except 
P  are  now  known.  The  method  of  determining  P  is  ex- 

S  X  S  X 

plained  in  a  later  section;  refer  to  Eq.  (5.15).  Once  P 

S  X 

has  been  found,  Eq.  (4.13)  fixes  f(Mx)  whereupon  Eq.  (4.3) 
fixes  itself.  The  solution  of  Eq.  (4.3)  for  M  when 

f (M)  is  known  involves  an  iterative  procedure  which  is  ex¬ 
plained  later  in  this  section. 

From  Eqs .  (4.6),  (4.7),  (4.8)  and  (4.10)  we  also  find  that 

(4.14) 


P  f  (M  ) 
sx  x 


=  P  f(M  )  =  P  ,  f(M.) 
sy  y  s3  3 


Assuming  P  and  M  known,  this  relation  along  with  Eq.  (4.3) 

S  X  X 

fixes  P  when  M  is  specified;  it  also  fixes  f(M.J  and 

sy  y  3 

when  Pg3  is  specified.  The  methods  of  calculating 
and  Pg2  are  explained  later;  refer  to  Eqs.  (6.8)  and  (7.14). 

Fig.  4.1  is  a  rough  sketch  of  the  function  f(M)  as  defined 
by  Eq.  (4.3)  .  This  function  vanishes  at  M  =  0  and  at  M  =  °°  .  • 
The  maximum  value  occurs  at  M  =  1  and  equals 

(Y  +  1) 

(4.15) 


max 


2(y  -  1) 


.Y  +  1 

It  is  clear  from  the  figure  that  in  the  range  0<f<  f  , 

max 

Eq.  (4.3)  has  two  positive  real  roots,  one  subsonic  and  the 
other  supersonic.  Usually  the  supersonic  root  must  be  discarded 
because  it  is  not  consistent  with  a  further  constraint  imposed 
by  the  second  law  of  thermodynamics. 
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In  order  to  solve  Eq.  (4.3)  for  M  when  values  of  f 

and  y  are  specified,  we  utilize  an  iterative  procedure  based 

on  Newton's  method.  Let  M  denote  the  nth  trial  value  of 

n 

M  .  We  then  define  a  function  F(M  )  and  its  derivative 

n 

t 

F  (M  )  as  follows, 
n 


F(Mn)  =  f 


1  +  M  ' 

2  n 


(Y  +  1) 
I  2  (y  -  1) 


-  M 


F*  (Mn>  *  Mn  1  + 


(3  -  y) 

2Ty  TT 


n 


-  1 


(4.16) 


(4.17) 


The  (n  +  l)at  trial  value  of  M  may  now  be  taken  as 


F(M  ) 

M(n  +  1)  =  Mn  -  -^JL 


F  (Mn) 


(4.18) 


The  cycle  of  calculations  defined  by  Eqs.  (4.16),  (4.17) 
and  (4.18)  is  repeated  until  the  result  for  converges  to 
a  stable  value  at  the  desired  level  of  accuracy. 

Thus 


M  =  £im  M 

n 

n  -*•  00 


(4.19) 


It  is  readily  apparent  that  the  value  of  M  so  obtained 
satisfies  Eq.  (4.3)  for  the  prescribed  values  of  f  and  y  . 
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5.  Entropy 


Consider  the  case  shown  in  Fig.  2.1(a).  The  specific 
entropy  sx  of  the  discharged  gas  may  be  expressed  in  the 
form 


where  Px^  and  P  ~  are  the  partial  pressures  of  the  two  com¬ 
ponents  of  the  gas  mixture.  Symbols  PQ  and  Tq  denote  the 
pressure  and  temperature  of  the  ambient  atmosphere.  The  en¬ 
tropy  of  each  component  is  assigned  the  value  zero  at  this 

reference  state  P„  ,  T 

o  o 

Incidentally,  the  primary  fluid,  being  condensible,  might 
well  exist  only  in  the  liquid  state  at  the  reference  condition 
PQ  ,  Tq  .  Nevertheless,  so  long  as  it  is  in  the  gaseous  state 
at  station  x  ,  Eq.  (5.1)  may  still  be  used. 

The  partial  pressures  of  the  two  components  are  propor¬ 
tional  to  the  respective  mol  fractions.  Hence 


(5.2) 


(5.3) 
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Upon  substituting  Eqs.  (5.2)  and  (5.3)  into  (5.1)  and 


simplifying  we  may  obtain  the  result  in  the  form 


where 


(5.4) 


=  R3  £n  R3  -  x  Rx  An  (x  R^)  -  (  1  -  x)R2  An  [(1  -  x)R2J 

(5.5) 


It  can  be  shown  that  this  quantity  sm  represents  the 
entropy  increase  associated  with  the  mass  mixing  of  the  two 
different  gases. 

Moreover,  since  the  entropy  of  the  stagnation  state  s^ 

is  by  definition  identical  to  that  of  the  corresponding  static 

state  x  ,  we  may  replace  T  and  P  in  Eq.  (5.4)  by  T 

XX  sx 

and  Psx  ,  respectively.  Thus  we  obtain 


sx 


sm  +  Cp3 


(5.6) 


In  previous  studies  of  this  series,  we  have  considered 
only  the  special  case  in  which  the  same  fluid  is  employed  for 
both  the  primary  and  secondary  streams.  In  that  case  we  have 


(5.7) 


C.=C,=C,  =  C 
pi  p2  p3  p 


(5.8) 


YX  =  Y2  =  y3  =  Y 


(5.9) 


Moreover,  in  this  special  case  there  is  no  mass  mixing,  so 
that  we  must  set 
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0 


(5.10) 


LJiui  yguuu 


S 

m 


Consequently,  Eq. 

s  =  0  + 
x 


(5.6)  now  simplifies  to 


(5.11) 


Next  reverting  to  the  case  of  two  different  gases,  let 
At  denote  a  time  interval  during  which  unit  mass  crosses 
station  x  .  During  this  same  interval  x  units  of  mass  cross 
station  1  and  (1  -  x)  units  of  mass  cross  station  2.  The 
entropy  which  leaves  across  station  x  must  equal  the  entropy 
which  enters  across  stations  1  and  2  plus  the  entropy  increase 
caused  by  the  mass  mixing.  Thus 


(5.12) 


(5.13) 


(5.14) 


Upon  eliminating  sx  between  Eqs.  (5.6)  and  (5.12)  ,  we 

find  that  also  cancels  from  the  result.  It  is  then  a 

m 

simple  matter  to  solve  for  the  exit  stagnation  pressure  in  the 
form 


x  s. 


(1  -  X)s2 


(5.15) 
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Thus  the  sequence  defined  by  Eqs .  (5.13),  (5.14)  and 

(5.15)  now  fixes  P  .  Next  referring  back  to  Eq.  (4.13), 

sx 

we  can  calculate  the  value  of  f  (M^)  ;  finally,  from  Eqs. 

(4.16) ,  (4.17),  (4.18)  we  can  calculate  itself.  Once 

P  '  T  *  T  ,  and  Mv  are  known,  it  is  a  simple  matter  to 

SX  SX  S  w  X 

calculate  the  corresponding  conditions  Px  ,  from  the 

usual  perfect  gas  relations. 
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6 .  Momentum 

Consider  the  idealized  one  dimensional  constant  area 
ejector  shown  in  Fig.  2.1(b).  The  momentum  equation  for  this 
system  may  be  written 

(P1  +  P1  V12>A1  +  (p2  +  p2  V22)A2  =  (py  +  py  Vy2)Ay  (6,1) 

The  mass  conservation  relation  for  this  system  can  be 
expressed  in  the  form 


Pi  A,  ^ 


P2  A2  V2 

“(T  ' -“T 


=  p  A  V 

y  y  y 


(6.2) 


Let  us  divide  the  three  terms  of  Eq.  (6.1)  by  the  corres¬ 
ponding  three  terms  of  Eq.  (6.2).  Notice  that  the  areas  can¬ 
cel  out  of  the  result.  In  this  way  we  obtain 


<P1  ♦  P2  V> 

pT7! 


(P, 


+  (1  -  X) 


+  P2  V22) 
p2  v2 


(P  +  p  V  2) 

y  y  y 

p  v 

y  y 


(6.3) 

Temporarily  dropping  the  station  subscript,  we  next  deve¬ 
lop  the  typical  term  of  Eq.  (6.3)  in  generalized  format  as 
follows . 


(P  +  PV2) 
PV 


(1  +  VM2) 

M  /T  /T 
s' 


(6.4) 


Observe  that  the  pressure  P  cancels  from  the  result. 
It  is  convenient  to  define  the  auxiliary  function 
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(6.10) 


ssBraKHesag 


Fig.  6.1  discloses  the  existence  of  a  supersonic  root  over 
the  range  gw<g<gmax.  This  root  should  not  be  accepted  until 
a  check  calculation  is  made  to  determine  whether  the  result  is 
consistent  with  the  second  law  of  thermodynamics.  The  subsonic 
root  is  found  by  retaining  only  the  negative  sign  before  the 
radical  in  Eq.  (6.8). 

Once  My  has  been  found  from  Eq.  (6.8),  the  corresponding 
value  of  Pgy  may  be  found  from  Eqs.  (4.3)  and  (4.14).  This 
suffices  to  fix  all  properties  at  states  y  and  sy  . 
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7.  Availability 

The  entropies  s1  and  s 2  of  the  two  input  streams 
were  defined  earlier  in  Eqs.  (5.13)  and  (5.14).  The  corres¬ 
ponding  steady  flow  availability  functions  with  respect  to  an 
ambient  atmosphere  at  pressure  PQ  and  temperature  Tq  may 
be  written 


h  -  cPi  (Tsi  -  V  -  To  si 


li<,  =  C  ,  (T  ,  -  T  )  -  T  s, 
2  p2  s2  o  o  1 


(7.1) 

(7.2) 


Consequently  the  total  available  energy  entering  the  sys¬ 
tem  becomes 


^m  ~  X  ^1  +  (1  “  x) ^2  (7.3) 

The  respective  entropies  of  the  streams  leaving  at  stations 
x  ,  y  and  3  are 


where  sm  has  been  previously  defined  in  Eq.  (5.5). 
The  corresponding  availabilities  of  the  streams 


(7.4) 


(7.5) 


(7.6) 


leaving 


at  stations  x,  y  and  3  may  therefore  be  written 


x  Cp3 

(Ts3  - 

V 

-  T 

o 

s 

X 

(7.7) 

y  “  cp3 

(Ts3  - 

To) 

-  To 

s 

y 

(7.8) 

3  =  Cp3 

(Ts3  - 

To} 

-  T 

o 

S  3 

(7.9) 

Now  consider  the  losses  of  availability  listed  below. 
By  substituting  the  above  expressions  for  the  quantities  on 
the  left  and  simplifying  we  readily  obtain  the  expressions 
shown  on  the  right.  Thus 


('1'  -  'PJ  =  T  s 

m  x  o  m 


(7.10) 


<¥ 


x 


R,  T 
3  o 


(<P 


x 


*3> 


R,  T 
3  o 


(7.11) 


(7.12) 


These  results  are  very  significant.  Eq.  (7.10)  defines 
the  loss  of  available  energy  caused  by  the  mass  mixing  of  the 
two  different  gases.  Eq.  (7.11)  defines  the  further  loss  of 
available  energy  caused  by  the  momentum  mixing  in  an  idealized 
ejector.  Eq.  (7.12)  defines  the  augmented  loss  of  available 
energy  caused  by  the  momentum  mixing  in  the  real  ejector. 

It  is  useful  to  postulate  an  empirical  relation  between 
these  last  two  losses  which  is  shown  below  on  the  left.  The 
expression  on  the  right  then  follows  from  Eqs.  (7.11)  and 
(7. 12)  .  Thus 


(♦ 

7? 


x 

x 


.  _  ln  <*./*.*> 
E "  u.ytt 


(7.13) 


.-.JZmkn 
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We  term  nE  the  ejector  effectiveness.  Its  value, 
which  is  always  less  than  unity,  must  be  estimated  from  appro¬ 
priate  test  data  on  ejectors. 

If  we  treat  nE  as  known,  Eq.  (7.13)  fixes  Pg.j  •  The 
solution  is  simply 

1_ 

We  can  next  find  ffM^)  from  Eq.  (4.14)  and  from 

Eqs.  (4.16),  (4.17)  and  (4.18).  This  suffices  to  fix  all  pro¬ 
perties  at  states  3  and  s3  thereby  completing  the  analysis  of 
the  ejector. 
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8.  Diffuser 


It  is  convenient  to  analyze  the  diffuser  in  terms  of  the 
concept  of  availability  in  steady  flow.  However,  we  wish  to 
deal  with  the  characteristics  of  the  diffuser  itself  and  these 
cannot  depend  on  the  arbitrary  values  Pq  ,  To  which  happen 
to  characterize  the  condition  of  the  ambient  atmosphere.  We 
can  achieve  our  objective  by  defining  the  availability  with 
respect  to  a  reference  state  P^  ,  which  characterizes 
the  diffuser  rather  than  with  respect  to  the  ambient  state 

Po  '  To  * 

Under  these  circumstances  we  may  write  the  appropriate 
availabilities  at  diffuser  inlet  and  outlet  as  follows 


(8.1) 


s4 


=  C 


P3 


<TS3  -  T3) 


(8.2) 


These  equations  make  use  of  the  fact  that  =  Ts3 

Also,  states  s3  and  3  are  at  the  same  entropy  so  that 


(ss3  -  s3>  *  0  *  Cp  -  *3 


(8.3) 


From  Eqs.  (8.1)  and  (8.3)  we  conclude  that 


*s3  =  Cp3  (Ts3  -  T3>  " 


(8.4) 
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a*' 


This  result  shows  that  the  available  energy  at  the 


diffuser  inlet  is  simply  the  inlet  kinetic  energy  itself. 

The  loss  of  availability  through  the  diffuser  may  be 
found  by  subtracting  Eq.  (8.2)  from  (8.1).  The  result  is 


s4 


W 


Aip 


r3  t3 


We  now  define  the  diffuser  effectiveness  as 


so  that 


'Ui  -  ih  A 

ISJ _ !s4 

v  ^s3  / 


At|»  _  n  _  » 
*s3 


h3  T3  In 
ep3  lIs3  -  T3> 


This  reduces  readily  to 


Solving  for  the  pressure  ratio  gives 


e 


(1  -  nD> 


y3  m32 

2 


(8.5) 


(8.6) 


(8.7) 


(8.8) 


(8.9) 


This  is  the  result  required.  It  fixes  Pgij  when  Pgg  , 
and  nD  are  specified. 


iMSES  - ^ 


W  W 
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Assuming  the  area  ratio  A3/A4  to  be  arbitrarily  spec¬ 
ified,  we  can  find  from  the  continuity  relation.  Refer¬ 

ring  to  Eq.  (4.8)  we  may  write 


itl  ^  ■“ 


Ps3  A3 
A3  Ts3^y3 


f(M3)  = 


Ps4  A4 
/*3  *sP^3 


f(M4)  =  *4 


(8.10) 


Consequently 


f(M4) 


(8.11) 


Eq.  (8.11)  fixes  f(M4)  .  Then  M4  follows  in  the  usual 
way  from  Eqs.  (4.16),  (4.17)  and  (4.18).  The  result  suffices 
to  fix  all  properties  at  states  4  and  s4  .  Thus 


(8.12) 


*  WWV 


(8.13) 


*3 


P4  "  Ps4/<Ps4''P4>  <8-15’ 

The  configurations  studied  in  earlier  work  did  not  incor¬ 
porate  a  diffuser.  For  the  purpose  of  comparing  the  results 
of  the  present  analysis  with  corresponding  earlier  work,  it  is 
desirable  to  be  able  to  eliminate  the  effects  of  the  diffuser 
in  certain  cases.  This  can  be  accomplished  in  the  present 
analysis  simply  by  setting 
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1 


(8.16) 


^  nj.  If.  WM 


nD  = 


and 


(8.17) 


It  is  easy  to  confirm  that  under  these  circumstances 


Ts4  -  ^sS 


(8.18) 


P  .  =  P  , 
s4  s  3 


(8.19) 


m4  =  m3 


A4  =  A3  etc. 


(8.20) 

(8.21) 


and  all  diffuser  effects  disappear. 
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9.  Electrical  Power  Section 

The  power  conversion  process  in  an  EHD  duct  is  usually 
treated  on  the  basis  either  of  constant  area  or  of  constant 
static  state.  Inasmuch  as  electrical  power  output  per  unit 
mass  is  small,  the  numerical  differences  between  the  results 
computed  by  these  two  methods  is  negligible.  For  definite¬ 
ness  in  this  analysis,  however,  we  assume  constant  area. 

The  electrical  power  that  can  be  obtained  from  an  EHD 
duct  of  constant  area,  negligible  change  of  density  and  optimum 
length  can  be  estimated  from  the  one  dimensional  version  of 
Poisson's  equation  which  governs  the  electrical  field.  The 
solution  is  well  known  and  will  not  be  derived  here;  a  detailed 
derivation  may  be  found  in  Ref.  (3) .  The  essential  result  may 
be  written  in  the  form 

pe  m  7  SEb  a4  v4  (9-1> 

where 

Pe  =  gross  electric  power  output,  watts 

e  =  permittivity  of  medium 
-12 

=  8.854  x  10  farad/m  (for  any  gas) 

E^  =  dielectric  strength  of  medium,  volts/m 
=  area  of  duct  (constant) 

=  velocity 
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According  to  the  test  data  analyzed  in  Ref.  4,  the  dielec¬ 


tric  strength  is  well  approximated  by  the  expression 

Eb  =  CB3R3P  4 


(9.2) 


where  CB3  Is  a  characteristic  constant  of  the  medium.  The 
data  show  that  Eq.  (9.2)  applies  to  air  or  steam  up  to  about 
10  atmospheres  pressure  with 


CB3  ±  9.49  x  103  m2  uK/cmb  (9.3) 

In  this  study  we  assume  that  Eq.  (9.2)  can  be  extrapolated 
up  to  about  100  atmospheres  pressure.  We  also  assume  that  an 
expression  of  the  same  form  applies  to  other  media  besides  air 
or  steam  but  that  each  medium  has  its  own  characteristic  value 
of  the  breakdown  constant  CQ. 

In  the  two-fluid  system,  the  separate  values  of  the 
breakdown  constants  CBl  and  Cb2  are  usually  known,  but  the 
breakdown  constant  Cb3  of  the  resulting  mixture  is  seldom 
known.  In  the  absence  of  adequate  test  data  bearing  on  this 
point,  we  tentatively  assume  that  Cb3  can  be  estimated  from 
the  hypothetical  relation 

CB3  =  l^iCfli  +  (!  “  x)R2Cb2^  (9.4) 

Upon  substituting  Eq.  (9.2)  into  (9.1)  and  dividing 
through  by  the  mass  flow  rate,  we  find  the  electrical  work  out¬ 
put  per  unit  mass  of  fluid  in  the  formi 


P 


e 


*4 


pi 


p4 


Ts5  ) (9.5) 
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By  use  of  standard  perfect  gas  relations,  Eq.  (9.5)  can 


be  changed  to  a  more  useful  form.  In  this  connection  set 


<Ts4  -  Ts5>  "  ATs 


(9.6) 


and 


where 


6  = 


£  CB3  Po 


(9.7) 


2 

PQ  =  ambient  pressure,  N/m 
Tq  =  ambient  temperature,  °K 


Notice  that  3  represents  a  dimensionless  version  of 
the  electrical  breakdown  constant. 

Eq.  (9.5)  now  yields  the  important  result 


The  corresponding  gross  electrical  work  output  per  unit 
mass  of  fluid  is  then  simply 


W*  =  Cp3  ATs  (9.9) 

Unfortunately,  the  known  value  of  e  and  the  typical 
experimentally  measured  values  of  Cg  are  very  small.  Con¬ 
sequently,  the  dimensionless  breakdown  constant  3  is  ex¬ 
tremely  small  compared  to  unity.  Hence  the  drop  in  stagnation 
temperature  AT  and  the  gross  electrical  work  per  unit  mass 

S 

W*e  are  also  typically  very  small.  This  is  a  very  basic  and 
serious  limitation  on  the  performance  that  can  be  obtained  from 
an  EHD  generator. 

Treating  the  process  4>5  as  isentropic,  we  may  write 
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(9.12) 


This  fixes  f(Mg)  whereupon  Mg  may  be  found  in  the 
usual  way  from  Eqs.  (4.16),  (4.17)  and  (4.18).  However,  con¬ 
vergence  should  now  be  very  rapid  since  we  may  take  as  a 
first  approximation 

«  M.  (9.13) 

5  4 

This  solution  now  suffices  to  fix  all  properties  at 
states  4  and  s4  . 
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10.  Condenser/Separator 

The  primary  fluid  which  leaves  the  condenser/separator 
at  station  6  is  assumed  to  be  compressed  liquid  at  known 
static  pressure  and  at  known  static  temperature 

T.  =  T.  .  The  kinetic  energy  of  the  liquid  at  station  6  is 
assumed  to  be  negligible. 

Temperature  Tg  must  be  equal  to  or  greater  than  the 
ambient  temperature  Tq  in  order  to  satisfy  the  requirements 
of  a  heat  balance  on  the  condenser/separator.  No  attempt  is 
made  in  this  analysis  to  formulate  this  actual  heat  balance; 
instead  Tg  is  simply  treated  as  a  given  or  known  quantity. 
However,  it  may  be  of  interest  later  to  investigate  the 
effects  of  varying  Tg  on  the  overall  performance  of  the  EHD 
system. 

The  secondary  fluid  which  leaves  the  condenser/separator 
and  enters  the  ejector  at  station  2  consists  primarily  of  a 
low  molecular  weight  noncondensing  gas.  It  also  contains  some 
secondary  vapor,  but  the  mass  fraction  of  this  vapor  is  treated 
as  negligible  in  this  analysis.  Thus  the  secondary  fluid  at 
station  2  is  treated  as  dry  gas.  The  kinetic  energy  at  station 
2  is  not  negligible  and  must  be  taken  into  account  in  the  ana¬ 
lysis  . 

In  analyzing  the  the rmo dynamic  process  in  the  condenser/ 
separator,  it  is  again  useful  to  employ  the  concept  of  avail¬ 
ability  in  steady  flow.  The  availability  function  is  always 
defined  with  respect  to  some  suitable  reference  state.  Nor- 
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T  of  the  ambient 
o 


mally  this  is  fixed  by  the  conditions  PQ  , 

atmosphere.  In  some  cases,  however,  we  seek  to  characterize 

the  performance  of  some  component  of  the  system  in  a  manner 

which  is  independent  of  the  ambient  conditions  P  ,  T  .  In 

o  o 

such  cases  it  is  useful  to  choose  as  a  reference  state  some 
condition  that  actually  occurs  in  the  component  in  question. 

We  used  such  a  procedure  to  good  advantage  in  analyzing  the 
diffuser  and  in  defining  the  diffuser  effectiveness  nD  •  We 
now  use  it  again  in  analyzing  the  condenser/separator  and  in 
defining  a  friction  loss  coefficient  cf  which  characterizes 
the  operation  of  this  component.  For  this  purpose  we  choose 
as  the  appropriate  reference  state  of  zero  availability  the 
static  conditions  Pg  ,  Tg  in  the  condenser/separator. 

It  also  clarifies  matters  further  if  we  choose  the  con¬ 
dition  Pg  ,  Tg  as  the  datum  state  of  zero  enthalpy  and  zero 
entropy.  In  this  way  we  ensure  that  our  analysis  in  this  sec¬ 
tion  makes  no  reference  whatever  to  the  state  P_  ,  T 

o  o 

Remember  also  that  P,  =  P-  =  P .  and  T,  -  T.  . 

6  2  1  6  2 

Thus  we  may  now  write 

s5  *  sm  +  Cp3  (ff)  -  *3  ln  (jf) 


(10.2) 


(10.3) 


The  quantity  sm  in  Eq.  (10.1)  is  the  same  as  previously 
defined  in  Eq.  (5.5).  It  expresses  the  entropy  increase  asso¬ 
ciated  with  mass  mixing. 
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The  availabilities  at  stations  5,  2  and  6  are  therefore 


*5 

'  S3 

<Ts5  - 

V  "  T6  S5 

(10.4) 

*2 

CM 

ii 

(Ts5  - 

o 

1 

* — » 

CM 

(10.5) 

*6 

=  h6  “ 

■  T6  s6 

=  0 

(10.6) 

The  irreversibility  I  of  the  process  which  occurs  in  the 
condenser/separator  with  respect  to  the  specified  reference 
state  Pg  ,  Tg  may  be  expressed  in  terms  of  the  above  avail¬ 
abilities  as  follows 


I  =  \p5  -  (1  -  x)  \l>2  ~ 


(10.7) 


We  postulate  that  this  irreversibility  can 
lated  with  the  inlet  kinetic  energy  in  terms  of 
friction  loss  coefficient  cf  as  follows 


I  =  c . 


_5_ 

2 


:f  Y3  R3  Ts5' 


M5 


also  be  corre- 
an  empirical 


(10.8) 


where  as  usual 


(10.9) 


Of  course  the  above  friction  loss  coefficient  cf  must 
be  estimated  from  suitable  experimental  data. 

Assuming  that  c^  is  a  specified  constant,  the  trial 
solution  obtained  for  an  arbitrary  trial  value  of  x  and 
arbitrary  fixed  values  of  the  various  other  input  parameters 
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will  not  necessarily  satisfy  Eq.  (10.7).  It  is  therefore 
convenient  to  represent  the  residual  error  in  this  equation 
in  the  form  below.  The  factor  Tg  is  inserted  on  the  left 
to  make  E(x)  dimensionless.  Thus 

R3  Tg  E (x)  =  *5  -  (1  -  x)  y2  -  0  -  I  (10.10) 


The  foregoing  relations  can  be  assembled  and  rearranged 
as  follows.  Let 


Then  the  final  energy  balance  equation  reduces  to 


E(x)  =  K  +  K.  x  (10.13) 

o  1 

By  holding  all  other  input  parameters  fixed  and  allowing 
x  to  take  on  a  range  of  possible  values,  it  will  normally  be 
possible  to  find  by  trial  a  value  of  x  such  that 

E(x)  =  0  (10.14) 

This  then  is  the  value  of  x  which  is  consistent  with 
the  other  specified  input  parameters  and  which  defines  a  pos¬ 
sible  EHD  thermodynamic  cycle. 
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Of  course,  if  the  various  fixed  input  parameters  are  not 

properly  chosen,  Eq.  (10.14)  cannot  be  satisfied  for  any  trial 

value  of  x  which  lies  in  the  useful  range  0<x<l  and  no 

solution  exists  under  these  conditions. 

If  a  solution  exists  it  can  be  found  by  some  suitable 

iteration  technique.  One  useful  method  is  to  set  E(x)  =  0 

in  Eq.  (10.13)  and  then  solve  this  equation  for  x  but  to 

relabel  this  as  x  .  Then 

new 


K 

X  =  - 

new  Kx 


(10.15) 


Assuming  that  conditions  have  been  so  specified  that  a 
solution  for  x  actually  exists,  and  assuming  that  an  initial 
trial  value  of  x  has  been  chosen  that  is  sufficiently  close 
to  the  true  root,  repeated  application  of  Eq.  (10.15)  should 
finally  yield  convergence  to  the  true  root  itself. 

Incidentally,  any  trial  solution  for  which  E(x)<0  re¬ 
presents  a  physical  impossibility  because  it  violates  the 
second  law  of  thermodynamics.  Any  trial  solution  for  which 
E(x)>0  satisfies  the  second  law  but  represents  a  condition 
in  which  the  losses  in  the  condenser/separator  exceed  the 
value  prescribed  by  Eq.  (10.8).  Only  the  case  for  which 
E(x)  =  0  satisfies  both  Eq.  (10.8)  and  the  second  law. 
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11.  Cycle  Efficiency 

Once  a  solution  has  been  found  such  that  E(x)  =  0  ,  it 
becomes  a  straightfoward  matter  to  calculate  the  correspond¬ 
ing  overall  cycle  efficiency. 

Because  of  the  pressure  drop  through  the  boiler/super¬ 
heater  ,  the  pressure  P7  =  Pg7  at  boiler  inlet  is  slightly 
higher  them  the  stagnation  pressure  Pgl  of  the  primary 
fluid  at  the  ejector  inlet.  We  assume  that  the  ratio  Pg^/Pg7 
is  a  specified  constant.  Hence  with  Pgl  specified,  Pg7 

is  also  known.  Also  the  pressures  Pg  =  P^  are  known. 

* 

Consider  the  ideal  pump  work  wp  done  per  unit  mass  of 
primary  fluid.  For  the  present  consider  the  hypothetical  case 
of  a  reversible  pump.  Also  note  that  density  changes  across 
the  pump  are  negligible.  Hence  we  may  write 


(Ps7  -  V 
p6 


(11.1) 


The  ideal  gross  electrical  work  output  per  unit  mass  of 

mixed  fluid  has  earlier  been  established  in  Eq.  (9.5).  It  is 

★ 

denoted  by  symbol  wq  . 

Let  us  now  denote  the  net  useful  electrical  work  output 

per  unit  mass  of  primary  fluid  by  symbol  w  t  •  It  may  be 
*  * 

related  to  wg  and  wp  in  the  following  way 


w 


net 


Joule/kg 


(11.2) 


where  nx  ,  the  excitation  efficiency,  allows  for  the  small 
electrical  power  expended  to  excite  the  system  and  where  n  * 

r 
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the  pump  efficiency,  allows  for  irreversibilities  in  the  pump. 

The  heat  input  in  the  boiler/superheater  per  unit  mass 
of  primary  fluid  may  be  written 

* 

lin.  *  (hsl  '  h6>  -  ^  (11-3> 

where  hg  is  the  enthalpy  of  the  primary  liquid  at  conditions 
Pg  ,  Tg  as  listed  in  suitable  tables  of  properties. 

The  overall  thermal  efficiency  of  the  cycle  can  now  be 
calculated  from  the  simple  formula 


n 


c 


w 


net 


‘in 


(11.4) 


Both  n  and  w  .  are  useful  parameters  which  charac- 
c  net 

terize  the  overall  thermodynamic  performance  of  the  EHD  sys¬ 
tem. 

Eq.  (11.2)  shows  clearly  how  the  ejector  serves  as  a 
kind  of  amplifier  which  increases  the  electrical  work  output 
per  unit  mass  of  primary  fluid.  This  effect  is  shown  by  the 
presence  of  the  parameter  x  in  the  denominator  of  the  first 
term.  Notice  the  beneficial  effect  of  a  low  value  of  x  on 
parameters  wnet  and  nc  .  Unfortunately,  the  value  of  x 
cannot  be  stipulated  independently  in  advance;  it  is  fixed 
by  the  other  specified  input  parameters  as  explained  in  the 
previous  section. 

* 

Notice  that  owing  to  the  fact  that  wg  is  normally 
very  small,  it  is  possible  in  some  circumstances  for  wnet 
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to  take  on  values  which  are  actually  negative.  When  this 
happens  it  means  that  electrical  power  output  is  insufficient 
to  drive  the  pump. 

The  relationships  developed  in  this  and  the  preceding 
sections  make  it  possible  to  carry  out  systematic  parametric 
studies  of  various  one-fluid  and  two-fluid  EHD  cycles,  with 
realistic  allowances  for  the  various  losses  that  occur.  Such 
studies  can  establish  optimum  design  parameters  and  perfor¬ 
mance  limits  under  various  circumstances.  These  results  in 
turn  can  finally  permit  informed  conclusions  to  be  drawn  con¬ 
cerning  the  ultimate  feasbility  of  this  general  type  of  EHD 
power  generator. 


12.  Calculation  Sequence  for  Two-Fluid  System 

In  this  section  we  summarize  the  various  equations  pre¬ 
viously  derived  in  the  approximate  order  in  which  they  would 
be  used  in  the  calculation  of  system  performance.  Also  listed 
are  the  initial  input  parameters  whose  values  must  be  specified 
in  order  to  start  the  calculation  and  various  further  input 
parameters  whose  values  must  be  specified  in  order  to  proceed 
with  various  successive  stages  of  the  calculation. 

Initial  Input  Data:  p]_  =  p2  =  P6  '  Ti  '  Psl  '  Tsl  '  (hsl  “  hl}  ' 

R  =  8315  Joule/kg°K 


*1 


*n(T  ,/T .)\ 

-  OTP 


(12.1) 


(h 


Tt 


si 


si 


(12.2) 


wi  -  (12-4> 

Further  Input  Data:  Wy  ,  WQ  =  W2  ,  T2  =  Tg  ,  pv(Tg) 


CD  = 


W 

lP2  -  W  J 


Verify  that  cd<<1 


(12.5) 
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“  A-1/A3)  /^s2\  1 1^3  ^*s3 

f(M*)  =  f  U  -  *>“  ra  VR2  Ts2  y3  f(V 


( Y  3  +  1) 


max  \  y  -i  + 


^3  ' 


Verify  that  f(Mx)  =  f<fmax 

For  first  approximation  set  Mx  =  M2 

Iterate  using  Eqs.  (12.31),  (12.32),  (12.33). 


V  - £  ( 1  +  Mn2) 


(y3  +  1) 

2  (y~  -“If 


O  -  Y,> 


F'(Mn)  =  f 


(y3  ♦  !)  (  ,[ Y3  -  1)  2V  ‘Y3  ‘  11 

■J-T-  Mrl? ‘  V  M„y 


(12.27) 


(12.28) 


(12.29) 


(12.30) 


(12.31) 


(12.32) 


F  (M  ) 


L(n  +1)  Mn  F'  (Mn) 


(12.33) 


Iteration  converges  to  fix  M 


M1  \  Y1  "  1  2 

g(M  )  =  - L—  \  1  +  (-i_ - )  m. 

(1  +  Yx  M^)  V 


(12.34) 


g(M2)  = 


(1  +  Yo  M-  ) 


1  +  (- 


Yo  -  1 


(12.35) 


g(My) 


(1  -  2y3  g^)  -  y  1  -  2  (y3  +  1)  g" 
1  -  Y-,  (1  “  2y,  g2) 


(12.36) 


(12.37) 


1  \  (Y3  "  1} 

Eq.  (12.37)  applies  provided  that  g  y1  gm  =  —  \ - 2 -  (12.38) 

y  3  i 

If  g  =  g*,  the  solution  reduces  to 


My  " 


(Y,  -  1) 


Once  My  is  known  from  Eq.  (12.37)  or  (12.35) 


f(My) 


My  1  + 


l)  2 
-  My 


(y3  +  1) 


lY,  - 


Psy  Psx  f (My) 
Further  Input  Data: 


f(Mx) 


Ps3  =  Psx  i  p 


£(m3,  .  *«V 

Verify  that  f(M,)  =  f<f,. 


(12.39) 


(12.40) 


(12.41) 


(12.42) 


(12.43) 


(12.44) 


(12.45) 


For  a  first  approximation  take  M3  = 

Iterate  using  Eqs.  (12.31);  (12.32),  (12.33). 
Iteration  converges  to  fix  M3  . 


p3  *  *.3'ir.3'*3> 


Further  Input  Data:  nD 


Ps4  *  Ps3  e 


'( 


1  -  n 


f  (M  ) 

4 


f 


f  (M3) 


Verify  that  f(M4)  =  f<fmax 

For  a  first  approximation  take  M4  =  M3 

Iterate  using  Eqs.  (12.31),  (12.32),  (12.33). 
Iteration  converges  to  fix  M4  • 


(12.46) 

(12.47) 

(12.48) 

(12.49) 

(12.50) 

(12.51) 

(12.52) 

(12.53) 


(12.54) 
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^s4  (^3  “  l)  2 

xf  =  1  +  V 


T4  =  Ts4/(Ts4/T4) 


P  T  (y3  *  l! 

ps4  _  _s£  J 

P4  "  T4 


P4  "  Ps4/(Ps4/P4> 


(12.55) 


(12.56) 


(12.57) 


(12.58) 


Further  Input  Data:  e  =  8.854  x  10  z  farad/m,  Cg-^  ,  Cg 3 r 
3  2 

(C0  =  9.5  x  10  m  “K/cmb  for  air  or  steam) 

•D 


CB3  =  R3  R1  CB1  +  ^  “  x)  R2  CB2 ^ 


£  CB3  Po 


(12.59) 


(12.60) 


6  (^Mfe)2 11 * 


(1  -  —■) 
rs4 


T  .  V(y3  "  L) 

Ps5  -  Ps4 


(12.62) 


(12.63) 


Verify  that  f(Mg)  =  f<fmax 

For  a  first  approximation  take  Mg  =  M4 
Iterate  using  Eqs.  (12.31),  (12.32),  (12.33) 
Iteration  converges  to  fix  Mg  . 


(12.64) 

(12.65) 

(12.66) 


P5  ’  Ps5/(Ps5/P5l 
Further  Input  Data:  cf 

—  Rg  2.n  Rg  ~  x  R^  £n  (x  R^) 


-  (1  -  x)  R2  £n  [(1  -  X)  R2] 


K  =  Jn 


s5 


+  (Cp3  ~  Cp2)  /  Js5 


-  1 


m 


Cf  Y; 


M. 


K,  -  SEifjsS  - 

1  R3  \  T6 


(12.67) 

(12.68) 

(12.69) 

(12.70) 

(12.71) 


(12.72) 

(12.73) 
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Verify  that  KQ<0 


(12.74) 


(12.75) 


Iterate  from  Eq.  (12.14)  until  *  converges. 
Further  Input  Data:  '  p6  '  (hsl  ’  hs' 

PS7  -  Psl/lFSl/Ps7' 

»;  -  v 


(12.76) 

(12.77) 


(12.78) 


RESULT  (12.79) 


(12.80) 


w 


net 


RESULT  (12.81) 


13.  OPTIMUM  VELOCITY  RATIO 

The  calculation  sequence  summarized  in  the  preceding  sec¬ 
tion  reveals  that  a  comple  cycle  calculation  requires  us  to 
specify  numerical  values  for  about  thirty  input  parameters. 

One  of  these  is  the  velocity  ratio  (V^/^)  =  y  .  The  calcula¬ 
tion  procedure  outlined  above  then  permits  us  to  find  the  cor¬ 
responding  value  of  the  mass  flow  ratio  x  ,  assuming  that 
there  exists  a  value  of  x  which  satisfies  the  governing 
thermodynamic  relations  for  the  specified  values  of  the  other 
input  parameters. 

If  we  now  change  the  value  of  y  but  hold  all  other  in¬ 
put  parameters  fixed,  we  can  repeat  the  above  procedure  and 
find  a  corresponding  new  value  of  x  ,  if  such  a  value  exists. 

It  is  evident,  therefore,  that  x  becomes  some  definite 
function  of  y  as  long  as  we  remain  in  the  domain  where  a 
real  solution  exists.  Of  course  all  other  dependent  variables 
of  the  cycle  are  also  functions  of  y  .  In  particular,  the 
overall  cycle  efficiency  nc  is  some  definite  function  of  y  . 
Moreover,  for  prescribed  values  of  the  other  input  parameters, 
there  will  be  some  definite  value  of  y  ,  let  us  designate  it 
as  the  optimum  value  Y0pt  >  which  yields  the  greatest  value 
of  ti  that  is  possible  under  the  specified  input  conditions. 

Our  basic  purpose  therefore  is  to  determine  for  any  prescribed 
values  of  the  other  input  parameters,  the  value  of  yopt  and  the 
corresponding  value  of  r)c)rnax*  Of  course  the  values  of  xQpt 

and  other  dependent  variables  are  also  of  interest. 
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Once  a  working  computer  program  is  available  to  establish 


values  of  yQpt  /  nc^max  '  xcpt  311(2  so  on  for  Prescribed 
values  of  the  other  input  parameters,  we  can  then  proceed  to 
study  the  effects  of  changes  in  these  other  parameters  and  to 
search  systematically  for  such  values  of  the  most  important 
input  parameters  as  will  yield  the  best  overall  performance 
of  the  system. 
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14 .  CALCULATION  SEQUENCE  FOR  ONE-FLUID  SYSTEM 

The  calculation  sequence  summarized  in  Section  12  simpli¬ 
fies  for  the  special  case  where  the  same  condensible  fluid  is 
used  for  both  the  primary  and  secondary  streams.  This  special 
case  is  summarized  below. 


Initial  Input  Data;  =  V  ^  =  Pg  , 


t2  =  t6 


Assuming  pressure  P^  and  quality  z ^  are  specified,  we 
may  find  from  tables  of  properties  of  the  saturated  fluid 


hi  =  hg  -  a  -  z1)hfg 


(14.1) 


S1  “  sg 


-  (1  -  Zj_) 


(14.2) 


Then  with  s-^  fixed  and  Pgl  specified,  we  may  also  find 
Tgl  and  hgl  from  tables  of  properties  of  the  superheated 


vapor. 


In  addition  with  Tg  specified  somewhere  in  the  range 
TQ  Tg  <_  T-^  ,  and  with  Pg  =  Pj^  specified,  we  may  find  fi-  f  the 
enthalpy  of  the  liquid  at  state  6  from  the  tables  of  properties. 

The  primary  fluid  can  now. be  approximated  by  an  "equivalent 
perfect  gas"  whose  properties  are  as  follows 


Y  ,  L  .  tn(Tsi/Ti>) 

l  l"<Psl/Pl>' 
.  <h3l  ~  hl> 

P  -  <Tsl  -  Tl) 


(14.3) 

(14.4) 


(14.5) 
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The  secondary  fluid  is  represented  by  the  same  "equivalent 
perfect  gas"  and  is  therefore  assigned  the  same  values  of 
Y  ,  Cp  and  R  . 

Further  Input  Data:  (V  /V  )  =  y 

-i.  ^ 


(14.6) 


(14.7) 


Verify  that  M2  <  1 


(14.8) 


(14.9) 

(14.10) 


ff)  -  (Ts2/T2I 


y/(y  -  1) 


Ps2  *  P2(Ps2/P2> 

Further  Input  Data:  x  (trial  value) 
T 


s  3 


=  x  Tsl  +  (1  -  x)Ts2 


f(Mx)  =  Mx 
f  (m2)  =  m2 


1  - 


1  + 


r  (y  +  1)/2(y  -  1) 


r(Y  +  1)/2(y  -  1) 


(14.11) 

(14.12) 


(14.13) 


(14.14) 


(14.15) 


lpslj  V  Ts2  rap-  (1  -  X) 


(14.16) 
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© 


VA2 


(l  +  al/a2) 

Further  Input  Data:  PQ  ,  TQ  (ambient  conditions) 

cs 


“Gfl  •  i 


sx 


f(‘V 


pc e 


po 

'  To 

( V 

1  Tc 

j  -  RJln 

(h 

1  To 

j  -  RJln 

p 

R?) 

Jin  (P 

sx^o^ 

XS^“  (l-x)  S2 


-  11  -  \fls3  f(M) 

sx  W  2 


1  A  . 
(1  -  x)  [] 


\(  Y  +  1)/2(Y  ~  1) 


max  yr  +  ly 

Verify  that  f (Mx)  =  f  <  fmax 

For  a  first  approximation  set  Mx  ~  and  iterate 

Iterate  using  Eqs.  (14.26),  (14.27),  (14.28). 


I  Y  -  1  2\(y  +  1)/2(Y  '  1} 

U  -  f  ( 1  + 


F(M 


f'(mj  =  f(1-^-i)Mn(i  +  (]L-jH1)Mh 


Mn 


n 


(3  -  y)/2(y  -  1) 


-1 


F(Mn) 

M(n  +  1)=  Mn  "  F'  (Mn) 


(14.17) 

(14.18) 

(14.19) 

(14.20) 

(14.21) 

(14.22) 

(14.23) 

(14.24) 

(14.25) 

(14.26) 

(14.27) 

(14.28) 
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Iteration  converges  to  fix  M  . 

X 


g(M,) 


g(M,)  = 


(1  +  ymJ) 

M2 

(1  +  yM^) 


'i+CL^r-±)  m; 


(14.29) 


(14.30) 


g(M  )  -  g  )  .  v  \  T 


g  =  g 


x  WT 3l  ,  (1  -  x)  1  / Ts2 

M7r\jT33  M  Ts3 


M2  =  (1  ~  2yg2)  -  n/1  -  2(y  t  1)  g2 
y  1  -  Yd  ~  2yg2) 

It  can  also  be  shown  that  in  the  special  case  where 
1  /Y  -.'"I 

m  -  —  v-i— 2 —  ,  the  above  solution  degenerates  to 


(14.31) 


(14.32) 


“y  = 


Once  My  is  found  from  Eq.  (14.32)  or  (14.33) 

v  -  1  21  "  {Y  +  1)/2(Y  "  1) 

f(My)  =  My  1  +  (X-^-A)M2 

. 

f  (II  ) 

P  =  P  * 

sy  sx  f (My) 


Further  Input  Data:  nv 


(14.33) 

(14.34) 

(14.35) 


Ps3 

Psx 

fe) 

f  (M3) 

-fel 

f  (Mx 

(14.36) 


(14.37) 


Verify  tnac  f(Il3)  =  f  <  (14.38) 

For  a  first  approximation,  take  ~  (14.39) 

Iterate  using  Eqs .  (14.26)  ,  (14.27)  ,  (14.28)  . 

Iteration  converges  to  fix  • 


—  =  i  i 


T3  =  Ts3/('rs3/T3) 


Y/(Y  “  1) 

pf  *  <WT3> 


P3  “  Ps3/(Ps3/P3) 

A3 

Furtner  Input  Data:  n..  >  - — 

™  Lf  iri  • 


?S4  -  Ps3  e 


-a  -  n„i  (y  t ii/2) 


£(M4»=  £  f«"3> 


Verify  tnac  f(M. )  =  f  <  f 

J  4  nax 

For  a  first  approximation  take 


*  =  (3 


M- 


Iterate  using  Eqs.  (14.26),  (14.27),  (14.28) 
Iteration  coverges  to  fix  >1^  . 


Ts4  "  Ts3 


■34 


1  +  nj 

2  4 


(14.40) 

(14.41) 

(14.42) 

(14.43) 

(14.44) 

(14.45) 

(14.46) 

(14.47) 

(14.48) 

(14.49) 
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■  WW'V 


(14.50) 


»P 

*4 


Ps4 

P4 

-  (T34/T4)Y/(Y  ■ 

1) 

(14.51) 

P4 

=  PS4'/(PS4/P4> 

(14.52) 

Note : 

If  it  be  desired 

to  omit 

the  diffuser, 

merely  set 

nD  =  1  and 

A3/A4  =  1  .  The 

result 

is 

Ts4 

—  rr> 

*s3 

(14.53) 

T4 

=  T3 

(14.54) 

Ps4 

=  Ps3 

(14.55) 

P4 

"  P3 

(14.56) 

M4 

=  M3  e  tc 

(14.57) 

and  all  effects  of  a  diffuser  disappear  from  the  result. 

-12 

Further  Input  Data;  e  =  8.354  x  10  farad/m  for  all  gases 

Cg=  9.5  x.  10^  m^  °K/cmb  for  air  or  steam 


n 

p 


AT 

s 


p 


si 


2 

eC:P0 
3  O 


=  6  h_r_i) 

2y 


=  Ts4  -  ATs 


"  Ps4(Ts5/T 


(14.58) 


v-l/(Y-l) 

Mjj  (14.59) 

(14.60) 

(14.61) 


fu3) 


(14.62) 
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Verify  that  f(M-)  =  f  <  fmax 
For  a  first  approximation  take  ~  . 

Iterate  using  Eqs.  (14.26),  (14.27),  (14.28) 
Iteration  converges  to  fix  M~  . 

/T__\  r.  ,\ 


(14.63) 

(14.64) 


1  tT~) 

-  1  +  M5 

(14.65) 

\  o 

T  - 

j 

-  "3i/(TSi/i5) 

(14.66) 

(W 

-  (Ts5/T5)Y/<Y  '  11 

(14.67) 

P3 

-  ps5/lps5'/t'5) 

(14.68) 

Further  Input  Data:  cf 

Note:  For  a  single  fluid  cycle  sm  =  0 


=  £n 


(ft  -iHhr-tft 


M 


K1 

E(x) 

xnew 


T  - 
si 


(Y  -  1) 


=  K  +  K,  x 
o  1 


K 


-  1 


5  (14.69) 

(14.70) 

(14.71) 

(14.72) 


Return  to  Eq.  (14.13)  and  iterate  until  x  converges  and 
E  (x)  =0  . 
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-  .iiAmU 


Further  Input  Data:  (Ps]/P37)  /  Pg  /  (hgl 
P 


s7 


=  Psl/(P8l/P87> 


h6)  '  nx  '  np 

(14.73) 


w’ 


=  C  AT 
P  s 


(14.74) 


w 


<Ps7  -  P6> 


(14.75) 


V/ 


n  w 
x  e 


net 


w* 

nP 


RESULT  (Joule/kg)  (14.76) 


^■in 


(hsl  "  V 


w* 

-  _£ 

np 


(14.77) 


-  Wfiet 
'  ^in 


RESULT 


(14.78) 
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15.  SAMPLE  CALCULATION  FOR  A  ONE  "FLUID  SYSTEM 

In  this  section  we  illustrate  numerically  how  the  calcu¬ 
lations  for  a  one  “fluid  system  may  be  carried  out  for  one 
arbitrary  but  typical  operating  point.  To  find  an  optimum 
operating  condition  would  require  the  repetition  of  similar 
calculations  for  a  large  number  of  such  trial  operating  points. 
No  attempt  is  made  in  this  section  to  find  such  an  optimum. 

We  arbitrarily  choose  to  calculate  a  case  which  is  in 
some  respects  similar  to  that  identified  as  case  B  in  the 


appendix  of  Ref.  (1) .  Thus 

Medium 

,  we  have 

for  both  of  these 

Steam 

Maximum  pressure 

Psl 

1500  psia 

Quality  at  state  1 

Z1 

0.95 

Ejector  effectiveness 

nE 

0.9 

Friction  factor 

Cf 

0.1 

Dielectric  strength 

CB 

9,490  m2  °K/cmb 

Diffuser  characteristics 

No  diffuser 

Excitation  efficiency 

nx 

1 

Pump  efficiency 

np 

0.9 

Pressure  ratio  across 
boiler 

Psl/ps7 

1 

Nevertheless,  despite  these  similarities,  these  two  cases 
are  not  identical.  Case  B  treats  the  densities  at  states  1,  2, 
x,  y,  3,  4,  5  as  equal.  The  present  analysis  takes  the  true 
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density  changes  at  these  states  fully  into  account.  Also  it 
is  expedient  in  the  present  compressible  analysis  to  define 


the  parameters  nE  and  c^  somewhat  differently  than  they 
are  defined  in  Ref  (1). 

Ref.  (1)  also  shows  that  a  comparatively  good  performance  is 
obtained  by  choosing 

pi  =  p2  =  p6  =  600  psia  (15.1) 

and 

V, 

TT  =  y  =  4.25  (15.2) 

v2 

Therefore  we  arbitrarily  choose  the  same  values  for  the 
present  numerical  example,  despite  the  fact  that  it  is  not 
entirely  equivalent  to  case  B  of  Ref  (1) .  We  also  choose  to 
set  T6  *  T1  (15.3) 


The  above  conditions  now  suffice  to  fix  the  following 

properties.  These  are  obtained  from  Keenan  and  Kayes  steam 

tables  and  are  listed  on  page  10A  of  Ref.  (1). 

T  .  =  1115.9  °R 

s  1 

Tx  =  T2  =  Tg  =  946.2  °R  (15.3) 


(hsl  -  h^)  =  78.9  Btu/lbm 

(hgl  -  hg)  =  673.9  Btu/lbm 

We  designate  ambient  conditions  as  follows 

pQ  =  14.696  psia  =  1.013  x  105  H/m2  (15.4) 

TQ  =  520  °R  =288.9  °K 
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The  properties  of  the  "equivalent  perfect  gas"  may  now  be 


computed  as  follows 


Y  =  1  - 


*n<Tsl/Tl) 


si/cr 


•  (- 


An (1115 . 9/946 . 
An (1500/600) 


2196 


(15.5) 


C_  = 


<hsl  *  hl> 


78.9 


p  (I  -  T1)  (1115.9  -  946.2) 


=  0.4649  Btu/lbm 


(15.6) 


T>  _  (Y  “  1)  p  _  0.2196  a  r  a  rt  \ 
R - Y -  Cp  1.2196  (0<4649) 


=  0.08371  Btu/lbm  "R  =  65.15  ft  Ibf/lbm  °R 


(15.7) 


The  following  constants  will  prove  useful  for  later  appli¬ 


cation 


V  ~  1  =  0.1098 


(Y  -  1) 


=  4.5537 


^-t1  =  1.1098 


=  5.0537 


(15.8) 


5-5537 


13_Z_ Xl  =  4  0537 

2(y-l) 


Also  the  dielectric  breakdown  constant  can  be  expressed  in 
dimensionless  form  as 


ECVo 


(8.854  x  10~12)  (9.49  x  103)  (1.013  x  10*5) 


(288.9) 


=  0.9678(10) 


(15.9) 


The  Mach  numbers  at  the  ejector  inlet  are 


M1  " 


2  ,Tsl  _ 
(Y  -  1)  T,  1} 


2  ,1115.9 

0.2196  l94.6i  ~  X) 


=  1.2781 


(15.10) 


M2  = 


M 


if  M1  -  I7T5  d-2781) 


=  0.3007 


(15.11) 


The  following  functions  are  calculated  for  later  use 

v  .  ,  o  -  (y  +  D/2(y  -  l) 

f (M1)  =  Mx(l  +  Y  i~  Mi  ) 

=  1.2781(1  +  0.1098  (1.2781)  2)  "  5,50537 

=  0.5552  (15.12) 


f(M2) 


M2(l  + 


Y  -  1  M  2.- (y  +  1)/2(Y  -  1) 
2  2  ' 


7  -5  0537 

=  0.3007(1  +•  0.1098(0.3007)  )  ' 

=  0.2861 


f 

max 


9  (Y  +  1)/2(y  -  1)  9  5.0537 

<y4 T>  - 

0.5907 


g(M1) 


MJi  +  M? 

(1  +  Y  Mx2) 

1.2781  Vl  +  0. 1098  (1. 2781)2 
(1  +  1.2196  (1.2781)  2) 


(15.13) 


(15.14) 


0.4639 


(15.15) 


=  0.3007  l/l  +  0.1098  (0.3007)  2 

(1  +  1. 2196(0. 3007)2) 

=  0.2722  (15.16) 

Stagnation  conditions  at  station  2  are  calculated  next. 

(Ts2/Ti)  =  1  +  (Y  ^  ■1)M22  =  1  +  0.1098  (0.3007) 2 

=  1.0099  (15.17) 

T2  =  T2(Ts2/T2)  =  946.2(1.0099) 

=  955.6  °R  (15.18) 

(Ps2/P2}  =  (Ts2/T2)Y/(Y  1)  =  (1 • 0099)  5  *  5537 

=  1.0562  (15.19) 

Ps2  =  Ps(Ps2/P2)  =  600-°  d-0562* 

=  633.7  psia  (15.20) 

In  this  particular  case  the  entropies  at  stations  1  and  2 
happen  to  be  equal.  They  can  be  calculated  either  from  the  stag¬ 
nation  temperatures  and  pressures  or  from  the  static  temperatures 
and  pressures.  The  latter  choice  gives: 

S1  =  s2  =  c^n^/T^  -  R  In  (p1/PQ) 

=  0.4649  Jln(946. 2/520.0)  -  0.08371  Jin  (600 . 0/14 . 696) 

=  -0.0322  Btu/lbm  °R  (15.21) 

Further  Data  Input:  x  =  0.472  (trial  value) 

The  above  value  of  x  was  found  by  previous  trial  calcula¬ 
tions  not  shown  here.  The  calculations  which  follow  verify  that 
this  value  of  x  does  indeed  satisfy  the  energy  balance  equation 
of  the  system  to  within  the  round-off  error. 
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Ts3  =  x  Tgl  +  (1  -  x)Tg2  =  0.472(1115.9)  +  0.528(955.6) 


=  1031.3  °R 


,Ps2.  .  Tsl  f(M2> 
psi  \*.2  W 


(15.21) 


_  0.472  ,633.7  ,  \  1115. 9 ,0.2861,  _  „ 

-  Q ;'Sg8'  ^ i 5 o 6 . 0'  '•  Qgy- z~ ( n  etcv  -  0.2103  (15.22) 


Al/A2  _  0.2103 
(1  +  A1/A2)  1.2103 


=  0.1738 


(15.23) 


Since  s^  =  s2  in  this  case,  Eq.  (14.20)  simplifies  to 


P  T 

*n(p££)  =  |[Cp  4n(^)  -  sL] 


=  0,-gaT/I  [0*4649  ( 526 /&'3)  +  °*0322]  =  4.1876  (15.24) 

p,Y  *  p„04*1876  =  14 . 696e  4,1876  =  967.9  psia  (15.25) 

sx  OS 


f(Mx)  = 


(i  -  ax/a3)  ps2  |ts3 


=  (1  o.528726)  'Ufc*'  \ fer  (0-2861)  "  °'3045  <15,26> 


Solving  for  Mx  from  Eqs.  (4.16),  (4.17),  (4.18)  gives 


Mx  =  0.3225 


(15.27) 


For  an  ideal  ejector 


si  «  (1  -  x)  ,  Ts2 


9(V-  +  MtJ 

-  (»  iMS)'1  -  0-3413  (!*.*.» 


2  =  (1  -  2yq2)  -  1/  1  -  2(y  +  1) g2 

'  1  -  y(1  -  2yg2) 
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[1  -  2.4392(0  .3418) *] 


2]  -  (l  -  4. 


4392 (0.3418) 


=  0.1658 


1  -  1.2196  [1  -  2.4392  (0 . 3418) 2 ] 
My  =  0.4071 


M 


f(M  ) 

y 


[i  + 


TTT  2  (Y  +  D/2(y  -  1) 
1 - -  M„2] 


2  y 
0.4071 


[1  +  0.1098(0.1658) ]^-0537 


=  0.3716 


(15.29) 


(15.30) 


f(Mx)  0.3045 

Psy  ~  Psx  f  (My)  =  967,9  (  0".'37 16)  =  793,1  Psia 


(15.31) 


p„..  l/n. 


'  s3 


,793.1. 


p-’  -  PSX(P^I  °  *  967-9  (5f777> 


1/0.9 


S3 


f  (M3)  =  (p^)f(Mx)  =  (0.3045)  =  0.3799 


=  775.7  psia  (15.32) 

(15.33) 


Solving  for  M3  from  Eqs.  (4.16),  (4.17),  (4.18)  gives 


l3 

=  0 

.4183 

In 

the 

absence  of  a  diffi 

s4 

= 

Ts3  =  1031.3  °R 

s4 

= 

Pg3  =  775.7  psia 

4 

= 

M3  =  0.4183 

4} 

= 

f(M3)  =  0.3799 

(15.34) 


(15.35) 


The  drop  in  stagnation  temperature  through  the  working 


section  is  given  by  the  following  expressions. 


(15.36) 


0.9678(10) 


,0.1098,  ,775.7  ,  ,  520.0, 

'nrrw  'uiw  (mo’ 


,  4.5537 

[1  +  0.1098  (0.4183) 


2 


-3 

1.0722(10) 


,  (15.37) 

ATs  =  TS4(ATS/TS4)  =  1031.3 (1.0722) (10)  =  1.1057  °R 


Notice  that  this  temperature  drop  is  only  about  1.1  °R  or 
about  one  tenth  of  one  percent  of  the  absolute  stagnation  tempera¬ 
ture! 


T  _  =  T  .  -  AT  =  1031.3  -  1.1  -  1030.2  °R 

s5  s4  s 


Ps5  ■  Ps4<Ts5/Ts4> 


Y/<Y  -  1) 


_  tic  -?  ,1030.2, 

"  775*7  (lff31.3} 


5.5537 


=  771.1  psia 


(15.38) 


(15.39) 


(15.40' 


(15.43) 

i 


To  verify  the  energy  balance  in-  the  condenser/separator  we 


calculate 


Y 


K°  - 


’s5, 


CjrY  T 


4n(^2i)  -  — {jjr-)  Mj 


=  £n( 


771. ls 
600. 


5.5537  4n( 


1030.2, 
946.2  } 


0.12196 ,1010.5 
2  1 546 . 2“ 


)  (0.4210) 


2 


=  -  0.2330 

# 


(15.44) 


K1  =  Ty  -~3T  -  !)  -  5.5537  (94§°y  -  1)  =  0.4930  (15.45) 

The  final  energy  balance  now  gives 
E(x)  =  Kq  +  K1  x  =  -  0.2330  +  0.4930  (0.472) 

=  -  0.0003  (15.46) 


This  verifies  that  the  energy  balance  is  satisfied  to  within 
round-off  error.  It  also  confirms  that  the  value  x  =  0.472  is 
the  mass  flow  ratio  that  corresponds  to  the  various  independent 
input  conditions  stipulated  in  this  sample  calculation. 

Now  that  the  value  of  x  has  been  confirmed,  we  may  complete 
the  thermodynamic  cycle. 

Since  the  pressure  drop  across  the  boiler  is  taken  as  negli¬ 
gible  in  this  example 

Ps7  =  Psl  =  1500-°  Psia  (15.47) 

The  gross  electrical  work  output  per  unit  mass  is 

w  *  =  C  AT  =  0.4649(1.1057)  =  0.5140  Btu/lbm  (15.48) 

e  p  s 

The  ideal  pump  work  is 
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(Ps7  -  V 


144  (1500.0  -  600.0) 

778.26  49.75 


=  3.3472  Btu/lbm 


(15.49) 


The  net  useful  work  of  the  cycle  is  therefore 


n  w 
x  e 


(1)  (0.5140)  3.3472 

0.472  "  0.9 


=  -  2.6301  Btu/lbm 


(15.50) 


Notice  that  the  net  work  of  the  cycle  happens  to  be  negative 
at  this  particular  operating  point.  In  other  words  the  electrical 
power  output  is  insufficient  to  drive  the  pump  under  these  condi¬ 


tions  . 


The  net  heat  input  to  the  cycle  is 


V  -  id- W 


-  670.2  Btu/lbm 


(15.51) 


Finally,  the  overall  thermal  efficiency  of  the  system  at  this 
operating  point  is 


net  _  -2.6301 


=  -  0.003924 


(15.52) 


Although  the  performance  at  this  arbitrary  operating  conditions 
is  very  poor,  the  fact  remains  that  there  exists  some  (as  yet 
unknown)  value  of  the  parameter  V]/V2  =  ¥  for  whioh  the  efficiency 
tv  reaches  its  maximum  possible  value. 

C 


This  optimum  value  of  y  can  be  determined  by  further 
systematic  trial  calculations  similar  to  the  above  sample 
calculation. 

The  foregoing  calculations  have  established  the  values  of 
T  ,  P  and  M  at  all  stations.  The  corresponding  values  of 
T  ,  P  ,  p  and  s  at  these  locations  can  then  be  determined 
from  the  following  relations. 


T 

=  1  +  (— — ^ — “) M2  =  1  +  0.1098M2  (15.53) 

T  =  T  / (T  /T)  °R  (15.54) 

s  s 


p„  T  Y/(Y  -  1)  5.5537 

=  (^)  =  (^> 


P  =  PS/(PS/P)  psia 

P  =  It  "  (^TT?}  T  lblR/ft3 

Ts  ps 

s  =  C  ln  (=2.)  -  RlnigS-) 


T  Pc 

=  0.4649&n(g-^  'q)  -  0.08371Jln(r|--^g-)  Btu/lbm  #R 


(15.55) 

(15.56) 

(15.57) 

(15.58) 

(15.59) 


The  calculation  of  the  above  properties  has  been  completed 
in  this  way  and  the  results  are  summarized  in  Table  15.1.  These 
same  states  are  shown  on  a  temperature  entropy  diagram  in  Fig.  15.1. 

Notice  from  the  tabulation  that  in  passing  through  the 
ejector,  the  fluid  density  p  changes  from  1.402  lbm/ft3  to 

3 

1.524  lbm/ft  ,  an  increase  of  about  8.5%  .  In  passing  through 
the  working  section,  the  density  changes  from  1.524  lbm/ft3  to 
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TABLE  15.1  THERMODYNAMIC  PROPERTIES 


4649  Btu/lbm 


TEMPERATURE- ENTROPY  DIAGRAM  OF  THERMODYNAMIC  STATES 


3 

1.515  lbm/ft  ,  a  decrease  of  well  under  1%  .  These  data  support 
the  principle  that  it  is  permissible  to  neglect  the  density  changes 
across  these  components  for  purposes  of  simplifying  the  analysis. 
This  is  particularly  justified  because  at  present  the  uncertainity 
in  the  value  of  the  empirical  ejector  effectiveness  nE  far 
exceeds  the  small  error  introduced  by  neglecting  these  minor  density 
variations.  Moreover,  the  present  sample  calculation  applies  to  a 
velocity  ratio  4.25  which  is  an  extreme  case.  For 

efficient  ejector  operation  the  ratio  V1/V2  should  be  only 
slightly  greater  than  unity  and  under  these  circumstances  the  den¬ 
sity  change  across  the  ejector  is  expected  to  be  far  less  than 
8%  .  Thus  there  is  little  substance  in  the  argument  that  an 
analysis  which  neglects  these  minor  compressibility  effects  must 
be  seriously  in  error.  Such  an  argument  runs  counter  to  the 
fundamental  principles  of  fluid  mechanics.  A  more  correct  appraisal 
of  the  situation  is  that  the  inclusion  of  these  detailed  compressi¬ 
bility  effects  can  be  expected  to  improve  the  accuracy  of  the 
analysis  to  some  degree,  but  that  they  should  not  be  expected  to 
change  the  overall  results  in  any  radical  way. 
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16.  Nomenclature 


B 


E(x)  - 

F<Mn>- 
f(M)  - 

fmax 
g(M)  - 


’max 

9co 

h 

I 

Ko 

K1 

M 

ih 

P„ 


qin 

R 

R 


area 

speed  of  sound 

characteristic  breakdown  constant  (Eqn  9.3) 

specific  heat  at  constant  pressure 

friction  loss  coefficient  -  condenser 

breakdown  electric  field  strength 

Eqn.  10.13 

Eqn.  4.16 

Eqn .  4.3 

Eqn.  4.15 

Eqn .  6.5 

Eqn .  6.7 

Eqn ,  6.9 

Eqn.  6.10 

enthalpy 

irreversibility 

Eqn.  10.11 

Eqn.  10.12 

Mach  number 

mass  flow  rate 

gross  electrical  power 

pressure 

heat  input 

gas  constant 

universal  gas  constant 
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s 

s 

T 

V 

W 

W€ 

WE 

w. 

X 

X 

y 

z 


m 


net 


new 


entropy 

entropy  of  mixing 

temperature 

velocity 

molecular  weight 

gross  work  out  per  unit  mass 

ideal  pump  work  per  unit  mass 

useful  electrical  output  per  unit  mass 

m1/m3 

Eqn.  10.15 

VV3 

quality  of  vapor 


B 


* 

<u 


Eqn.  9.7,  dimensionless  breakdown  constant 
ratio  of  specific  heats 
change  of 

permittivity  of  the  medium 

overall  cycle  efficiency 

Eqn.  7.13,  ejector  effectiveness 

Eqn.  8.6,  diffuser  effectiveness 

pump  efficiency 

excitation  efficiency 

gas  density 

steady  flow  availability 

mass  ratio  of  condensible  primary  to  dry  secondary 
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Subscripts 

0  -  ambient  conditions 

1  -  primary,  entrance  to  ejector  (exit  from  boiler/superheater) 

2  -  secondary,  entrance  to  ejector  (exit  from  separator) 

3  -  mixture,  exit  from  ejector  entrance  to  diffuser 

4  -  mixture,  exit  from  diffuser  entrance  to  generator 

5  -  mixture,  exit  from  generator  entrance  to  condenser 

6  -  primary,  exit  from  condenser  entrance  to  pump 

7  -  primary,  exit  from  pump  entrance  to  boiler 

f  -  saturated  liquid 

G  -  gas  or  secondary 
g  -  saturated  vapor 
n  -  next  value 
s  -  stagnation 
v  -  vapor 

x  -  exit  station  of  hypothetical  minimum  loss  device 
y  -  exit  station  of  ideal  ejector 
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